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ABSTRACT 

In this paper we study 21-cni absorption spectra and the corresponding emis¬ 
sion spectra toward bright continuum sources in the test region (326°< I < 333°) 
of the Southern Galactic Plane Survey. This survey combines the high resolution 
of the Australia Telescope Gompact Array with the full brightness temperature 
information of the Parkes single dish telescope. In particular, we focus on the 
abundance and temperature of the cool atomic clouds in the inner galaxy. The 
resulting mean opacity of the HI, < k >, is measured as a function of Galactic 
radius; it increases going in from the solar circle, to a peak in the molecular ring 
of about four times its local value. This suggests that the cool phase is more 
abundant there, and colder, than it is locally. 

The distribution of cool phase temperatures is derived in three different ways. 

The naive, “spin temperature” technique overestimates the cloud temperatures, 
as expected. Using two alternative approaches we get good agreement on a 
histogram of the cloud temperatures, TccoO corrected for blending with warm 
phase gas. The median temperature is about 65 K, but there is a long tail 
reaching down to temperatures below 20 K. Glouds with temperatures below 40 
K are common, though not as common as warmer clouds (40 to 100 K). 
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Using these results we discuss two related quantities, the peak brightness 
temperature seen in emission surveys, and the incidence of clouds seen in HI 
self-absorption. Both phenomena match what would be expected based on our 
measurements of < k > and Tcooi- 

Subject headings: 21-cm line, interstellar medium, radio astronomy 


1. Background 

Absorption by the 21 cm line of HI was first detected 48 years ago (Hagen, Lilley, 
and McClain, 1955), just four years after the line was detected in emission from interstellar 
hydrogen in the Galactic plane. As spectral line interferometers were developed in the 1960’s 
(Clark, Radhakrishnan, and Wilson 1963, Radhakrishnan et al. 1972a) it became possible to 
measure absorption spectra toward many more continuum sources at high and low latitudes. 
The striking result from these pioneering studies was that the Hi absorption spectra look very 
different from emission spectra in the same directions (Clark 1965, Hughes, Thompson, and 
Colvin, 1971, Radhakrishnan et al. 1972b). This can be understood as a result of the huge 
range of temperatures in the interstellar hydrogen (Field 1958, Field, Goldsmith, and Habing 
1969). Interstellar atomic hydrogen can exist in stable thermodynamic equilibrium either at 
warm temperatures, ~5000 to ~10,000 K, or at cool temperatures, <50 to ~150 K; these two 
phases are reflected in two regimes of excitation temperature for the 21-cm line. Generally 
this “spin temperature” is close to the kinetic temperature of the gas, but see Liszt (2001) for 
a detailed discussion of the exceptions to this rule. There is good evidence that interstellar 
hydrogen is present at temperatures intermediate between those of the classic warm neutral 
medium (WNM) and cool neutral medium (CNM) (Heiles 2001), but generally the different 
morphology of the emission and absorption spectra can be explained by a mixture of gas 
along the line of sight at widely different excitation temperatures (see Kulkarni and Heiles 
1988 for a deeper discussion of the phases of the interstellar medium). 

This paradigm for understanding the different structures of 21-cm emission and absorp¬ 
tion spectra as due entirely to variations in the excitation temperature has been challenged 
occasionally on the basis of the different physical volumes sampled by the two kinds of spec¬ 
tra. The problem is not so much the necessity of moving the telescope on and off the back¬ 
ground continuum source to obtain the absorption and emission spectra, respectively. The 
real problem is the different solid angles contributing. Emission at 21-cm is so widespread 
that the emission spectrum always corresponds to the beam area of the telescope, while 
the absorption spectrum comes only from the gas in front of the continuum source. The 
solid angles of the two can differ by several orders of magnitude, which raises the question 
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of whether the differences between emission and absorption could be due simply to small 
scale structure in the interstellar medium (Faison and Goss 2001, Deshpande 2000). To get 
around this problem requires an emission-absorption study using the smallest possible beam 
size for the emission and using background sources which are relatively large, if possible the 
same size as the beam. The new mosaic surveys of the Galactic plane, combining single dish 
and interferometer data, allow this for the first time. The Southern Galactic Plane Survey 
(SGPS, McGlure-Grifhths et ah 2000) and the Ganadian Galactic Plane Survey (GGPS, 
Taylor et ah 2002, Strasser et al. 2002) provide good quality emission-absorption spec¬ 
trum pairs using extended continuum sources at low latitudes, with a beam small enough to 
mitigate the effects of variations in the emission and absorption over small angles. 

In this paper we concentrate on a small, test region of the SGPS (McGlure-Grifhths et 
al. 2001) which has been mapped at relatively high resolution (FWHM 90"). Over the next 
year maps of the rest of the ~210 square degrees of the SGPS with similar resolution will 
become available, which will provide some 30 times the number of background sources as 
those considered here. The purpose of this paper is to test methods of analysing the emission- 
absorption spectrum pairs, and their interpretation. The next section (2) discusses the best 
way to obtain the spectra from the data. We then consider what the absorption spectra alone 
tell us about the opacity of the ISM in the inner galaxy (section 3). Then comes the tricky 
question of how best to combine the information from the emission and absorption spectra 
to estimate the spin temperatures of the cool clouds (section 4). We consider several htting 
techniques that parallel approaches used in past studies at intermediate latitudes. Finally 
we discuss the implications of the spin temperature distribution for observable quantities 
like the peak brightness temperature of the HI emission and HI self-absorption (section 5). 


2. How to Get the Emission and Absorption Spectra from the Data 

The HI 21-cm line is one of the only transitions in the entire electromagnetic spectrum 
for which emission and absorption are both relatively easy to detect from the same region. 
The two fundamental spectra are the brightness temperature of the emission, Tb{v), and the 
optical depth, r(n). Where there is no background continuum we observe Tb directly, and 
toward an extremely strong background source we observe r with negligable contribution 
from the emission, but toward most continuum sources we see a mixture of emission and 
absorption. This section deals with the optimum method for obtaining Tb and r from the 
survey data. 

There are always more weak background sources than strong ones; to get as many 
absorption spectra as possible in a given area we need to hnd the most effective way to extract 



the absorption spectra from the data, since the fainter continuum sources give absorption 
spectra of marginal quality. For this experiment, as for most studies of Galactic 21-cm 
emission and absorption, the limiting factor which sets the noise level in the absorption 
spectrum is not the radiometer noise but the precision with which we can subtract the 
emission in the directions of the background continuum sources, leaving the absorption only. 
This is a severe problem at low latitudes; the difficulty of doing accurate emission subtraction 
has made single dish low latitude emission-absorption studies very difficult, even for Arecibo 
(Kuchar and Bania, 1991). 

Interferometer studies offer such high resolution that the emission can always be elim¬ 
inated, at least toward the most compact background sources, but lacking the zero-spacing 
(total power) information, interferometer data alone cannot be used to hnd the emission 
spectrum. So for emission-absorption studies we need both single dish and interferometer 
data. The standard method for many years has been to measure the absorption spectrum 
with an interferometer and compare with an emission spectrum from a single dish telescope 
(Radhakrishnan et al. 1972c, Dickey et ah 1983, Mebold et al. 1982). In the SGPS we have 
carefully combined single dish (Parkes) and interferometer (ATGA^) data to achieve com¬ 
plete uv coverage for all angular scales larger than our synthesized beam size (90"). Thus for 
absorption studies we have the best of both worlds, with sufficient resolution to remove the 
emission in the direction of the continuum sources to get accurate absorption spectra at low 
latitudes, while still preserving the total power information which gives us the corresponding 
emission spectra. 

To obtain the most accurate estimate for the absorption spectrum we still must remove 
the emission as completely as possible from the spectrum toward the continuum source. 
This can be done optimally by considering the uv distribution of the continuum flux vs. 
the emission line flux. Generally the continuum is more compact spatially than the line 
emission, which translates to higher flux at longer uv spacings. The 21-cm emission follows 
a power-law distribution of flux vs. + ^ so that hltering in uv space allows us to 

reduce the emission in a predictable way. The hltering is performed by Fourier transforming 
the cube to uv space, then simply setting to zero all Fourier components within a given 
radius of the origin, then transforming back. This is illustrated on hgure 1, which shows 
the rms measured in the line channels vs. the uv radius of the hlter for a test region 1.42 
“square with little continuum emission (region 1 of Dickey et ah, 2001). The uv pixel size 
in this case is 1.4 (degrees)”^, meaning that a hlter radius of one uv pixel (which takes 
out just the zero-spacing point) leaves spatial structure with sizes smaller than ~0.7°, while 


^The Australia Telescope National Facility is funded by the Commonwealth of Australia for operation as 
a National Facility managed by the Commonwealth Scientific and Industrial Research Organisation 
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Fig. 1.— The rms of the emission after high-pass spatial hltering. The abscissa is the 
high-pass hlter cutoff radius on the uv plane, and the ordinate is the rms of the remaining 
emission. The 21-cm line emission is shown by the circles, the continuum is indicated by the 
stars. The continuum is made up of a combination of diffuse emission and discrete sources, 
while the line emission is mostly diffuse. Thus hltering has a stronger effect on the line 
channels than on the continuum. The pixel size in this case is 1.4 (deg)“^; the corresponding 
spatial wavelengths are shown on the top scale. 




a filter radius of ten pixels leaves only emission in structures smaller than ~4'. Figure 1 
shows that the rms level of the emission remaining after hltering is about 10 K for hlter 
width of one pixel, decreasing to 1.7 K for the ten pixel hlter. Also shown on hgure 1 is 
the continuum, using a region 1.42 “square centered on (l,b) = (331.43,0.56) that includes 
both diffuse continuum and a few bright sources. The continuum is only weakly affected 
by the hltering for hlter widths less than about hve pixels, but progressively more strongly 
attenuated for hlter widths ten pixels or more. On the basis of hgure 1 we choose a ten pixel 
hlter [passing spatial frequencies higher than (8.5 arc min)“^ which preserves angular sizes 
smaller than about 4'], since hltering more heavily begins to reduce the continuum as much 
as the line emission. Already with this hltering the total continuum hux of our background 
sources is less than it is on the unhltered map, since they are extended objects, thus we are 
beginning to lose signal to noise on the absorption spectra. We tabulate below these hltered 
continuum peak values since they determine the noise level in the absorption spectra, but 
these should not be used to estimate the hux densities of the continuum sources. 

In addition to removing the emission toward the continuum sources to obtain the ab¬ 
sorption spectra as described in the last paragraph, we need to interpolate the emission from 
nearby beam areas in order to estimate the emission spectrum which would be seen in the 
direction of the continuum source if there were no absorption. This we do by a spatial inter¬ 
polation on the unhltered cube, which includes information from all the short uv spacings. 
The technique is described by McClure-Grifhths et ah (2001); we construct for each spectral 
channel a bi-linear function htted to the pixels around the background source for which the 
continuum brightness is below a threshold set at 20% of its peak value on-source. We use 
the oh-source pixels so dehned both for the htting and to estimate the error of the htted 
function; the latter is given by the rms of the data minus the ht averaged over the set of 
oh-source pixels. That rms should be an overestimate of the likely error in the interpolated 
spectrum in the direction of the background source, since the area covered by the oh-source 
pixels is much larger than that of the source itself. The error envelope dehned this way is 
indicated above and below the emission prohles on hgures 2-8. The absorption spectra are 
constructed by averaging spectra toward the pixels for which the ratio of the continuum 
brightness to the continuum peak is above 80%, weighting by the square of this ratio. (The 
continuum map used to determine this weighting has also been spatially hltered in the same 
way as the spectral line cube.) Note that we do not need to subtract the interpolated emis¬ 
sion from the spectrum toward the continuum to obtain the absorption spectrum, since the 
spatial hltering process has already accomplished that step. This is the diherence between 
the analysis performed here and that of McClure-Grifhths et al. (2001). Comparing the 
spectra in hgure 8 of that paper with the corresponding spectra in hgures 2 and 3 below, 
we see that the ehects of emission huctuations have been attenuated by a factor of two to 
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three by the spatial hltering step described above. The accuracy of both the interpolated 
emission spectrum and the absorption spectrum is limited by this hltering and/or interpo¬ 
lation step, as neither can be observed directly without some pollution by the other. Thus 
the matched-hlter based on hgure 1 is of central importance. 

The brightest continuum sources in the SGPS test region are listed on table 1, which 
gives the source name, rms noise in optical depth {(Tt), integrated hux density in Jy, peak 
brightness temperature in K, and recombination line velocity, in columns 1-5. Columns 
6-10 give the values of various spectral integrals dehned below. All of these continuum 
sources are HII regions except for G328.42-1-0.22 which is a supernova remnant. The HII 
regions all have recombination lines detected in the survey of Caswell and Haynes (1987); 
the measured recombination line center velocity is an indication of the kinematic distance. 
In a few cases the HI absorption continues to velocities beyond (i.e. more negative than) 
the recombination line velocity. This may be an indication that the HII region is at the far 
distance corresponding to its velocity (G328.82-0.08 is an example). In some cases however 
the HI absorption spectrum shows a single deep component extending 5 to 10 km s“^ beyond 
the recombination line velocity, and then no absorption at higher velocities (G326.45-1-0.90 
and G326.65-1-0.59 are examples). These suggest local departures from the smooth circular 
rotation described by the rotation curve, perhaps due to the effect of the Sgr-Gar arm at 
about -50 km s“^. The recombination line velocities themselves may be offset from the HII 
region velocity due to asymetric expansion. 

In practice, peak continuum brightness temperatures of Tc > 50 K or higher are needed 
to give high quality absorption spectra (ut <0.03), but we include on table 1 all 13 sources 
with peak brightness above 20 K for the sake of the line integrals discussed in the next section. 
For the spin temperature discussion in the subsequent section we will restrict the analysis 
to the five sources with > 60 K. Gaussian fitting is done on the nine with > 30K. 
Spectra toward all the continuum sources listed on table 1 are shown on hgures 2-8, with 
the interpolated emission above and the absorption below. The intrinsic properties of most 
of these sources are discussed in detail by McGlure-Grifhths et al. (2001). 


3. Line Integrals and K{r) 

The absorption spectra toward the weaker continuum sources have so much noise due to 
emission fluctuations that we cannot compare the emission and absorption spectra channel- 
for-channel to derive spin temperatures for each cloud or line component, as we do for the 
spectra toward the stronger sources in section 4 below. But we can work with the velocity 
integrals even for the noisier absorption spectra, since emission fluctuations produce positive 
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Fig. 2.— Emission absorption spectrum pairs toward G331.35+1.07 (top) and G326.45+0.90 
(bottom). For both, the emission spectrum is at the top and the absorption spectrum 
is below. Both spectra have error envelopes, which are computed from the errors of the 
interpolation of the emission in the direction of the continuum source. The numbered boxes 
in between show the velocity ranges of the different blended lines. The rms noise in the 
absorption, coming primarily from the emission fluctuations, is indicated as cr^-, the peak 
continuum brightness temperature is indicated as Tc^max (in K), and the equivalent width, 
EW is indicated (in km s“^). Two values of the emission integral over negative velocities are 
shown, converted to units of 10^*^ cm“^, as discussed in the text. Also shown are Gaussian 
hts to the absorption (equation 9) and their residuals (data minus fit) for spectra toward 
continuum sources with Tc^max > 30 K. 
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Fig. 3.— Emission-absorption spectrum pairs for G326.65-1-0.59 and G327.99-0.09, presented 
as in figure 2. The absorption spectrum toward G326.65-1-0.59 was observed by Goss et al. 
(1972) in their pioneering study of Galactic 21-cm absorption. Our spectra are consistent 
with theirs within the uncertainties caused by emission fluctuations in their data. 
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and negative “psendo-absorption” (Radhakrishnan et al. 1972c) with no bias in the mean. 
The velocity integral of the optical depth spectrum is the equivalent width, EW, 


EW = [ t{v) dv (1) 

Jav 

where the range of integration, AV, is adjusted to cover a specihc interval along the line 
of sight using the rotation curve. Here we use the rotation curve of Fich, Blitz, and Stark, 
(1989) : 


_ 1.00746ReHe 

n + 0.017112^0 ^ ^ 

where H is the angular velocity of rotation around the Galactic center, Rgai is the Galacto- 
centric radius, and Gq and Rq are the corresponding values for the LSR at the solar circle 
(220 km s“^ and 8.5 kpc). 

Kinematic distances have a near-far ambiguity symmetric about the tangent point, but 
Rgai is unique, so we can assign to each radial interval in the inner galaxy a corresponding 
velocity interval, AV, and compute the corresponding EW. Dividing each EW by the line 
of sight length. As, through the corresponding annulus gives the mean opacity, < k >, by 


(k) 


EW 

As 



( 3 ) 


where n is the HI density and Tgp the excitation temperature of the A21-cm line, and the 
line of sight average denoted by the brackets, ( ), is taken over As which corresponds to the 
range of integration in velocity, AV, used to compute EW. Ch is the familiar conversion 
constant 1.823 x 10^^ ^ ^_i which converts to 5.9 x JAi , giving k units of km s“^ 

kpc“^. In these units k typically has values of a few to a few tens in the Milky Way disk. 
Note that < k > is the frequency integral of the more familiar (e.g. Spitzer 1977 eq. 3-2 
and 3-14), which has units (length)”^, although the range of integration may be only one 
spectral line channel. 
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Fig. 8.— Emission-absorption spectrum pair for G331.52-0.07, presented as in figure 2. As 
for G326.65-1-0.59 on figure 3, this absorption spectrum was first observed by Goss et al. 
(1972), with similar results to within their errors. 
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Results for (k) obtained by dividing the inner galaxy into radial bins of width ^-R© 
are shown on hgure 9. The histogram shows an arithmetic average over the various lines 
of sight through each annulus. Also shown on figure 9 as small circles are the individual 
channel values of r in the absorption spectra, converted to k by multiplying by the velocity 
gradient, given by the rotation curve (see Burton, 1988, for a review of the significance 
of the velocity gradient). 

The curves on hgure 9 indicate observational selection based on our noise level in r for the 
case of G328.42+0.22, one of our brighter sources. The lower curve shows the opacity which 
would result from optical depth equal to one sigma, as dehned by the emission huctuations 
discussed above. Weaker absorption than this is not detectable, so multiplying by the velocity 
gradient gives the corresponding lower limit for detectable opacity. Spectra toward fainter 
sources that give higher cTt will have higher lower limits. Points near or below this curve on 
hgure 9 are upper limits on k. On the high side, when the absorption lines are very deep, the 
noise prevents us from distinguishing between optical depths greater than Tmax = — In (2ar). 
The line saturates at this point (typically r ~ 3), and we set the optical depth to this value 
if e~^ < 2ar- Multiplied by the velocity gradient, this maximum detectable r gives an upper 
limit on the measurable k, shown for the case of G328.42+0.22 by the upper curve. Again, 
for spectra with higher noise level the upper limit is lower, as evidenced by chains of points 
at various levels corresponding to saturated lines. Points near the upper curve on hgure 9 
are thus lower limits on k. 

The velocity integrals of the absorption spectra are given on table 1, columns 6 and 
7. Golumn 6 gives the integral over the entire negative velocity range corresponding to the 
inner galaxy. Golumn 7 gives the integral over the restricted velocity range corresponding to 
the points on hgure 9, i.e. from zero km s“^ to the recombination line velocity. Whether the 
HII region is at the near or far distance, we can be conhdent that at least over this velocity 
range only the near-side gas can contribute to the absorption. For the Hi I regions we stop 
3 km s“^ short of the recombination line velocity, to avoid the deep absorption usually seen 
just beyond the HII region’s velocity. The SNR G328.42-1-0.22, is known to be beyond the 
solar circle on the far side of the Galaxy (Gaensler, Dickel, and Green, 2000); in this case we 
carry the velocity integration to the terminal velocity at this longitude. (Varying the 3 km 
s“^ ohset to 10 or even 15 km s“^ has minimal ehect on the opacities on hgure 9, although 
it reduces the number of lines of sight contributing to some of the annuli.) 

Figure 10 shows the geometry of the Galactic plane in the fourth quadrant (from 
McGlure-Griffiths et al. 2001). The sun’s location is assumed to be at {x,y) = (0,8.5), 
and we plot lines of sight at longitude 326°and 333°, the boundaries of the test region stud¬ 
ied in this paper. Annuli are plotted with radii from 0.4 to 1.0 R©, and the major spiral 
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Fig. 9.— Measurements of the 21-cm line opacity as a function of Galactic radius. The circles 
represent individual spectral channels, translated into optical depth per unit path length 
using the velocity gradient, as described in the text. The histogram shows the averages of 
many channels on many lines of sight for radial bins, weighting by line of sight length. The 
upper and lower curves indicate the observational limits due to saturation (top) and noise 
(bottom). 
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Fig. 10.— A face-on view of the Galaxy, with the longitude range of the SGPS test region 
indicated by two lines at longitudes 326°and 333°. The right panel has the radial velocities 
predicted from the rotation curve drawn as contours; the left panel shows some annular 
bins which are crossed by the lines of sight at the edges of the region. The annuli used for 
the averages shown on hgure 9 are much narrower than those shown here, and there are 
many more of them. Both hgures have the spiral arm model of Taylor and Gordes (1993) 
superposed. Gommon names for the spiral features are Norma (labelled “a”), Scutum-Grux 
(“b”), Sagitarius-Garina (“c”), and Perseus (“d”). 
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features from the model of Taylor and Cordes (1993) are reproduced. We see that the Norma 
spiral arm dominates the radial range 0.5 to 0.6 Rq near the tangent point, the Scutum-Crux 
arm at the near distance dominates radii 0.6 to 0.7 Rq, and Sagitarius-Carina dominates at 
about 0.85 Rq. There is a gap in the range 0.7 to 0.85 which is between arms. This appears 
to be reflected in a dip in k at about 0.8 Rq on hgure 9. Even in this dip, and at all other 
radii in the inner Galaxy, the opacity is larger than its solar circle value of about 5 km s“^ 
kpc“^. 

This narrow longitude range may not be representative of the entire inner galaxy, so 
we must wait for the full SGPS survey area to be studied before deciding conclusively that 
the opacity increases as we go to smaller Galactic radii. In the hrst quadrant, Garwood and 
Dickey (1989) found about the same result for the solar circle opacity, but smaller values 
for the inner galaxy. However, recently Kolpak et al. (2002) have made more extensive 
observations with the VLA in the first quadrant, and they hnd a trend similar to what we 
see on hgure 9. The difference may be due to the small numbers of lines of sight sampled in 
the Garwood and Dickey study, which may have preferentially sampled the regions between 
the spiral arms, which appear here to show opacities of 5 km s“^ kpc“^ or less, in contrast 
to the arm regions which can have < k >> 20 km s“^ kpc“^. 

If the increase in < k > by a factor of three to hve with decreasing Galactic radius from 
the solar circle to the molecular ring region at about 0.5 Rq is typical of the disk overall, we 
may ask whether it is due simply to an increase in the surface density of gas. The Hi surface 
density as traced by 21-cm emission surveys does not show an increase with decreasing Kgai 
in this longitude range (Burton 1988, £g. 7.15), but the molecular gas does. The strong 
variation of < k > tells us that the cool atomic gas, in contrast to the warm neutral medium 
that dominates the A21-cm emissivity, is much more abundant in the inner galaxy than at 
the solar circle. This underlines the role of the GNM clouds as an intermediate population 
between structures in the WNM and the much denser and colder molecular clouds. 

The integrals of the emission spectra are given on table 1, columns 8 and 9 in units of 
10^° cm“^. Golumn 8 gives the integral of the emission uncorrected for absorption, i.e. 

Nt = Ch J Tb{v) dv. (4) 

where we integrate over the full negative velocity range corresponding to the inner Galaxy. In 
column 9 we take the emission integral over the restricted velocity range AV corresponding 
to the absorption integral in column 7 (i.e. from zero to just short of the recombination 
line velocity). For this velocity range we can correct the values of emission for absorption 
to get a better estimate of the true column density assuming a one-phase medium (see next 
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section). That correction is given by 

Nc = Ch [ Tb{v) dv. (5) 

J Av i — e 

The effect of this correction is to increase the column density by a factor that is typically 
between one and two. Values of the ratio ^ for the same range of integration, 

AV are given on table 1, column 10. These correction factors show typical values of 1.4 
to 1.6, which is in agreement with similar correction factors found by Dickey and Benson 
(1982). Note that the Nf values on Table 1, column 8 are typically larger than on column 
9 because the range of integration is broader for the former than for the latter, the full inner 
Galaxy range for the former, and only to the recombination line velocity for the latter. 

Finally we can combine the velocity integrals of the emission and absorption to get a 
velocity averaged spin temperature. 


EW Ch 

which is the nominal value of the spin temperature needed to give the emission and absorption 
integrals if the gas were at a single temperature. This is given on table 1, column 11 for 
the velocity range corresponding to columns 7 and 9. Comparison of the emission and 
absorption spectra channel-by-channel shows that this value is far from realistic, even under 
the “one-phase” assumption. 


4. How to Combine the Emission and Absorption Spectra to find the Diffnse 

Cloud Temperatures 

Methods for combining the emission and absorption spectra to find the distribution of 
spin temperatures in the interstellar Hi can get complicated. The complications arise from 
blending of gas at different temperatures in the same spectral channel. Any line of sight at 
low Galactic latitudes will contain a mixture of several cool clouds which appear as Gaussian 
line components in the absorption spectrum, plus warm gas which is hardly visible in the 
absorption spectrum, but which contributes the bulk of the emission. The challenge is to 
combine the absorption and emission spectra in a way which separates these two thermal 
phases, giving an estimate of the temperature of the cool gas without the bias introduced 
by the warm gas emission. In this section we discuss and compare several mathematical 
techniques for this. 
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4.1. Method 1. The One-Phase Model 

The simplest way to combine the emission and absorption spectra is to ignore the 
blending and compute the excitation temperature for each velocity channel naively assuming 
all the gas at a given velocity is at the same temperature. This gives the simple formula: 


Tsp{v) = 


Tb{v) 
1 - 


(7) 


where Tb is the brightness temperature of the emission in velocity channel v, and t{v) is 
the optical depth of the absorption in the corresponding velocity channel. This gives the 
harmonic mean of the temperatures of the various neutral atomic regions which contribute 
to that velocity channel, i.e. 


T = 

sp 



( 8 ) 


where i is an index for the different regions on the line of sight contributing, each with 
column density Ni and excitation temperature Ti. The distribution of Tgp vs. number of 
channels for the nine sources with Tc >30K is given on hgure 11. We must assume that 
the temperatures resulting from equation 7 shown on hgure 11 are strongly biased to higher 
values than the true cool Hi temperature due to blending with the warm gas. 


4.2. Gaussian Fitting to the Absorption or to the Optical Depth 

Improved estimates of the cool phase temperature are based on the shapes of the features 
in the emission and absorption spectra. A simple approach is to ht the absorption spectrum 
with a number of Gaussian components. Low latitude A21-cm emission spectra do not 
generally look like the sum of distinct Gaussian components, but absorption spectra do, 
so this is a reasonable approach. When the Gaussian components blend together the best 
ht parameters become ambiguous, and for a complicated, low latitude spectrum even the 
number of Gaussian components needed is not clear. Thus there is no unique solution to the 
Gaussian htting problem. 

The htted parameters are obtained using an interactive program based on the Levenberg- 
Marquardt method of non-linear chi-squared minization (Press et ah, 1992). The numbered 
boxes above the spectra on hgures 2-8 indicate the velocity ranges dominated by the diherent 
line components (given on Table 2, columns 6 and 7). Weaker components may be hidden 
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Fig. 11.— Histogram of the temperatures derived from the naive, one phase model {Tgp, eq. 
7). The weighting is simply by number of channels, which roughly corresponds to weighting 
by volume rather than by mass. 
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below the stronger ones in these velocity ranges. Our detection limit for absorption features 
is thus a function of velocity, and is generally higher than that set by the noise. Note that 
these boxes do not represent velocity ranges over which the htting is performed. In cases 
where several line components are blended we perform the htting of all Gaussian parameters 
simultaneously over the full velocity range covered by all components, plus some baseline on 
either side. 

For deep absorption lines such as those found at low latitudes, it is preferable on physical 
grounds to ht the absorption spectrum with a sum of functions which are Gaussian in the 
optical depth, r, rather than in the observed absorption, e~'". We could do this by htting 
Gaussians to the log of the observed spectrum, but this amplihes the ehect of the noise 
for the deeper lines; for a deep line with noise we can get (1 — e~'^) > 1 which causes the log 
function to diverge. It is better to do the htting using for the model a sum of functions of 
the form 


(1 


g-rh) 


^ model 


1 


-E, 






( 9 ) 


where the parameters ro,i, vo,i, and ai for each line component i are varied to minimize 
the usual way (see table 2, columns 3-5). In spite of the double exponential this is no more 
difficult to ht to the data than a sum of ordinary Gaussians, and it provides a better match 
to the hat tops seen on the deeper lines. More importantly, the values of tq,* are signihcantly 
larger than the depths of the corresponding Gaussian hts, because of saturation. Thus the 
derived column densities for the cool gas are larger, often by a factor of two or more, than 
what would be obtained from simple Gaussian hts. The htted parameters are obtained using 
an interactive program based on the Levenberg-Marquardt method of non-linear chi-squared 
minization (Press et ah, 1992). The numbered boxes above the spectra on hgures 2-8 indicate 
the velocity ranges dominated by the diherent line components (given on Table 2, columns 6 
and 7). Weaker components may be hidden below the stronger ones in these velocity ranges. 
Our detection limit for absorption features is thus a function of velocity, and is generally 
higher than that set by the noise. Note that these boxes do not represent velocity ranges 
over which the htting is performed. In cases where several line components are blended we 
perform the htting of all Gaussian parameters simultaneously over the full velocity range 
covered by all components, plus some baseline on either side. 

We cannot estimate the cool phase temperature directly from the widths of the Gaus¬ 
sians, however they are htted, because turbulence and random motions on all scales have the 
ehect of widening the line signihcantly beyond the thermal width given by the Maxwellian 
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distribution of atomic velocities. So we need another approach to combine the emission and 
absorption spectra to get estimates for the temperatures of the cool phase gas and undo the 
blending of equation 8. The next simplest thing to do after the one-phase assumption of 
equation 7 is to assume the HI comes in two phases. 


4.3. Method 2. Two Phase Linear Least-Squares Fitting 

The two phase assumption attributes all the absorption to the cool clouds, while the 
emission spectrum includes the cloud emission plus emission from a warmer medium which 
is necessarily more broadly distributed in velocity. No assumption is made here about the 
spatial distribution of either phase, but the fact that optical depth spectra can be htted 
by a sum of Gaussians suggests that the cool gas is conhned to discrete structures. The 
assumption of two phases means that the observed brightness temperature is given by 

Tobs{v) = Tyjj + Twfi e + Tcooi (l — e Tcont e ^ (10) 

where the absorption comes from a cloud with temperature Tcooi and optical depth r, and 
the warm, optically thin gas is partly in front of the cloud, Tw,,/, and partly behind it, 

The last term provides for the presence of continuum emission, Tcont^ which is also behind 
the cloud, and hence is partially absorbed by it. Note that the contribution to the brightness 
temperature from the warm gas has nothing to do with its physical temperature; so long 
as it is optically thin its brightness temperature is proportional only to the column density, 
in the velocity channel width 5v^ as Ty^ = 

When we subtract the continuum from the emission spectrum to obtain the line bright¬ 
ness temperature, Tb{v), as a function of Doppler velocity, v, then the last two terms in 
equation 10 combine as 


Tb{v) — Tyjj + e + (Tcooi — Tcont) ~ ^ (H) 

This simplihes further if we consider the observed quantity (1 — e~'^) to be an independent 
variable that we rename x, so that 

Tb(v) = Tyjj T^^fe(l— x) -|- X (Tcooi — Tcont) (12) 

Next we make an assumption that is supported by observations at high and intermediate 
latitudes that the warm gas is broadly distributed in velocity relative to the widths of 
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the cool clouds. So over the velocity range of a single absorption line component we can 
approximate the warm gas contribution as a linear function of velocity, i.e. = Uq + aiv 
and Tu, 6 = 60 + biv or for the total 

Tw,f + = Co + Ciu. (13) 

with oo, Oi, bo, bi, cq, and Ci all constants. Then we have simply 


Tb{v) = (co + ciu) (1 - ex) + C 2 X (14) 

where C 2 is the cool cloud temperature reduced by the background continuum, C 2 = (Tcooi — 
Tcont) and e = ^— is the fraction of the warm gas which is behind the cloud (0 < e < 1 ). 
Here we will assume e is constant across the velocity width of each absorption line component, 
though in general it may be a function of velocity. Note that the continuum brightness, Tcont, 
is not the brightness temperature of the background source towards which we measure 
the absorption, but the diffuse continuum in the directions of the nearby pointings which 
give the interpolated emission spectrum. 

Given the observed absorption and emission spectra, we can perform a simple linear 
least-squares £t to the data in order to determine the values of Cq, Ci, and C 2 in equation 
14. The independent variables are v and x, and the dependent variable, Tb, we will call y 
for simplicity. The solution is obtained by minimizing the value of fhe same way as 
for polynomial fitting, which leads to a matrix inversion problem where the elements of the 
matrix are moments of various combinations of v, x, and y. For example. 


vxy 


1 

n 


Vi -Xi-pi 


(15) 


where the sum is taken over all spectral channels, i = 1 to n, covered by a distinct absorption 
feature (which presumably corresponds to a single cool phase temperature). The parameter e 
is carried through all the calculations explicitly, so that results can be found for any selected 
value of e. The advantage of this approach, compared with the Gaussian fits discussed above, 
is that there is no ambiguity in the results; the unique solution is found directly with no 
iteration and no need for a “first guess” input by hand. The equations needed to find this 
best fit are given in the appendix. 

This fitting method is illustrated in figure 12-13, which show plots of the emission 
brightness temperature, Tb against the observed absorption, 1 — e~'^. Different panels on 
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Fig. 12.— Detailed analysis of the variation of brightness temperature vs. optical depth. The 
top left panel shows at the top the emission-absorption spectrum pair toward G326.65-1-0.59 
as in hgure 3. The upper portions of the other panels show the same two spectra, but 
blowing up the velocity ranges indicated by the numbered boxes. The emission spectrum 
has superposed the least squares £t, and below the residuals (data minus fit). The linear 
function, T^,, is also shown with the emission spectrum. The graphs just below the absorption 
spectra are plots of Tb vs. x = (1 — e“^). The x’s on these plots show the individual channel 
values. The solid curve on the Tb vs. x plots shows the least squares htted function of 
the form of equation 14. The coefficients of this function are indicated for each absorption 
feature sub-plot by the equation at the bottom. 








3.— Plots of the spectral features toward G326.65+0.59, continued from figure 12 






figures 12-14 show separate line components toward G326.65-1-0.59, corresponding to the 
nnmbered boxes on the npper panel of hgure 3. The solid curves show the £t results, and 
the symbols show the measnred valnes in each velocity channel. In most cases the fits are 
qnite good, considering that there are only three free parameters, and typically 10 to 15 
independent pairs of measnred qnantities. 

The approach taken here is very similar to that of Mebold et ah (1997) and Dickey et al. 
(2000) in their stndies of the cool HI in the Magellanic Clonds. In those papers the fitting 
was done by hand by determining the “ridge-line” of each featnre measnred on the vs. 
r plane like those in figures 12 and 13. That approach works best for simple spectra withont 
blended absorption lines, as seen in the Magellanic Clouds and some high and intermediate 
latitnde Galactic directions. At low latitndes the many overlapping absorption featnres make 
it difficnlt to draw the ridge lines by hand, so the fitting techniqne described in this section is 
needed. The resnlts for nnblended lines are very similar to those derived nsing the ridgeline 
method. 

The fitted valnes of Cq, Ci, and C 2 tell us both the temperature of the cool cloud {Tcooi = 
C 2 + Tcont) and the linear £t to the warm gas emission across the velocity range of the 
absorption line component (eqnation 13). These fit resnlts are shown on table 2, columns 8 
throngh 10. For this analysis we use only the five spectra toward continnnm sonrces with 
Tc,max >75 K, for which the rms errors in the absorption are less than 0.025. From these fits 
we can derive new histograms of cool gas temperatnres. These are shown on figure 14, for 
three values of e=0.25, 0.50, and 0.75, which should bracket the true foreground - backgronnd 
distribntion. 

Figure 14 shows simply the number of spectral channels covered by clonds with a given 
valne of Tcooi- This is not a mass weighted average, nor a volnme average, bnt just the 
statistics of clouds sampled on the small nnmber of lines of sight toward the hve strongest 
continnnm sonrces. Since the sample is small, and gronps of channels come ont with the same 
clond temperature, the histograms have gaps at the lower temperatures, which are certainly 
unrealistic. These histograms must be considered merely representative of the distribntion 
of HI clond temperatnres, with the general conclusion that most clonds have cool phase 
temperatnres in the range 30 to 100 K, with occasional cooler valnes. 

The parameter e is nndetermined for any given clond. In general we may expect ronghly 
half of the warm-phase gas to be in front of the absorbing clond (at a given velocity) and 
half behind, which means e ~ 0.5 is a good gness. The lower panel, which shows the resnlt 
from setting e = 0.25 for all clonds, certainly nnderestimates the temperatnres of many 
clonds, since if only one qnarter of the warm phase emission comes from behind the clond, 
the cloud’s optical depth must be much higher, and hence its temperatnre cooler, then for 
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Fig. 14.— Histogram of Tcooi- The histograms show the number of spectral channels at each 
temperature, derived using method 2 (linear least squares). Results for three values of e are 
shown. The median values of Tcooi are 79 K, 65 K, and 51 K for e = 0.75, 0.5, and 0.25, 
respectively. 






30 



-100 -50 0 

Vlsr s ') 


Fig. 15.— Detail of the linear least squares fitting results for 0328.31+0.45. The upper 
spectrum (emission) shows both the data (Ts and the htted emission prohle, which are 
nearly indistinguishable, as well as the linear hts to the warm phase emission (cq + civ). 
These piecewise linear hts are not continuous, either with each other or with the emission 
spectra on either side of the absorption features. This motivates the third method (non-linear 
least squares) for hnding Tcoou which forces the warm phase models to be continuous. 
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e=0.5. Similarly, the upper panel, which assumes e = 0.75 for all clouds, must overestimate 
the temperature of most clouds. An exception is the cloud centered on -5 km s“^ toward 
G326.45+0.90. Fitting equation 14 to this results in a very cold temperature, 13 K assuming 
e = 0.75, just 2 K if e = 0.5, and -9 K for e = 0.25. (In the middle and lower panels of hgure 
14 this cloud drops off the left edge of the plot.) The latter two are clearly unphysical, so 
in this case we may conclude that e must be larger than ~ 0.6, and the cloud must be quite 
cold, strictly less than 25K, which is the value resulting from assuming e = 1. 


4.4. Method 3. Two Phase Non-Linear Least Sqnares Fitting 

The linear least-squares £t technique described by equations 10-15 and in the appendix 
does a good job of matching the observed values of Tb using the measured absorption values 
as the independent variable, but this technique ignores what we know about the warm gas 
emission at velocities on either side of the absorption line. Ideally the linear approximation 
to Ty, given by Cq-I-CiU should be continuous with the emission spectrum above and below the 
velocities covered by the absorption line, and in the case of multiple overlapping absorption 
lines, piecewise linear fits should be continuous between sequential velocity ranges. But this 
constraint is not incorporated into the fitting algebra, and typically it is not satisfied. An 
example is shown on figure 15, where the separate least squares fits toward G328.31-1-0.45 
are plotted together. The discontinuities between successive linear approximations to 
show that we are still missing some information that could improve the model of the HI 
temperatures. 

To force continuity in the model we have to return to a non-linear least-squares fitting 
technique, based again on the Levenberg-Marquardt method. This is a more complicated 
application, however, as we are now fitting a model to Tb{v) using two independent variables, 
V and a: = (1 — and a set of parameters. The parameters are C 2 = {Tcooi — Tcont) 

for each line component, and now instead of Cq and Ci to describe the warm gas emission 
for each velocity range, we have just one parameter for each contiguous box. For distinct 
velocity ranges (non-blended lines) the warm phase emission is fully determined by its values 
on either side of the absorption, and there is only one parameter to fit ( 02 ). For the case 
of overlapping lines, we can adjust the warm phase component arbitrarily at the boundary 
between the boxes, but in each box we stick to a linear interpolation for our estimate of 
the warm phase emission. Thus for blended absorption lines we are fitting both the cool 
phase temperature and the value of at each velocity edge where individual components 
overlap. The piecewise linear approximation to T^(u) is now fully specified by the values 
Tw{vi), where the set of uds are the boundaries of the velocity ranges represented by the 
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boxes on figure 15. This forces the piecewise linear model to be continuous from one 
component to the next. If we further require that be continuous with Tb on the first and 
last velocity edge, i.e. where the absorption drops below its detection threshold, then for a 
set of n overlapping line components we have to fit just n values of C 2 and (n — 1) values 
of Tw{vi), where i = 2, ...,n and Tu,{vi) and Tw{vn+i) are fixed. Figure 16 shows the same 
spectral pair as in figure 15 (G328.32+0.45), but now fitted using method 3. The estimate 
of the warm phase contribution is much better behaved. The residuals of the fit to Te^ are 
shown at the bottom. Given that there are fewer than two free parameters to fit for each 
absorption component, the agreement between the data and the fit is impressive. 

Figure 17 (lower panel) illustrates the distribution of cool phase temperatures resulting 
from this fitting for e=0.5, compared with a copy of the histogram for e=0.5 from figure 
14. The third method of fitting changes the results slightly, mainly by reducing the higher 
values of Tcooi and so narrowing the peak near 70K. There is still a tail reaching to low 
temperatures, and the median value of Tcooi is hardly changed (67 K vs. 65 K from method 
2 ). 

For a single, isolated absorption component, the warm gas level is determined by a 
linear interpolation between the values of Tb on the edges of the line, so there is only one 
parameter to fit, i.e. C 2 , that determines the cool gas temperature. This effectively reduces to 
the method of interpolation used by Kalberla, Schwarz, and Goss (1985) in their study of the 
emission and absorption in the vicinity of 3G147, and by Mebold et al. (1982) in their study 
of several high and intermediate latitude regions. If we were to assume that the absorption 
and warm phase emission components could all be fitted by blended Gaussians then we could 
derive the cool phase temperature directly. This is the approach taken by Heiles and Troland 
(2002a,b) in their Arecibo survey of absorption and emission at intermediate latitudes. In 
appendix 2 we make a detailed comparison of their results with the results of our Method 2, 
and find that the differences are minor. 


5. Discussion 

5.1. The Peak Brightness Temperature 

The mean opacity, < k >, derived in section 3 above plays a fundamental role in deter¬ 
mining the radiative transfer of the 21-cm line across the Galactic plane. At low latitudes, 
cool clouds with optical depths of one or more are common enough that some, but not all, of 
the velocity range in a typical direction will be covered by optically thick gas. How likely this 
is depends on a contest between < k> and the velocity gradient, along the line of sight. 




Fig. 16.— Detail of the non-linear least squares fitting results for G328.31-1-0.45, as in figure 
15. The piecewise linear hts for the warm phase emission are now much more reasonable. The 
data and the htted values for the total emission are also shown, but almost indistinguishable. 
The residuals (data minus £t) are shown at the bottom. 





10 20 50 100 200 

Tcooi (K) 


Fig. 17.— Comparison of the T^ooi distributions obtained with method 2 and method 3. 
Both histograms assume e=0.5. The upper histogram gives median T^ooi = 65 K, the lower 
histogram gives 67 K. Both are unweighted counts of spectral channels at each temperature, 
without regard to the total mass of gas in each channel. 





These two quantities have the same units (km s ^ kpc ^), so we can make a dimensionless 
ratio, p, dehned by 


K 
dv 
ds 

that depends on the local properties of the medium (the abundance of cool clouds) as well as 
on our particular vantage point (the longitude) and the shape of the rotation curve. In the 
extremely low velocity gradient case, as near longitudes 0°and 180°where ^ ~ 0, this ratio 
p is much greater than one, and the optical depth is high at all allowed velocities. In contrast, 
in the outer galaxy at a longitude where the velocity gradient is moderate, e.g. ±90 deg for 
which 1^1 ~ 10 km s“^ along most of the line of sight, high optical depths are rare, since 
p is less than one (< k > in the outer Galaxy is mostly below 5 km s“^ kpc“^). Since the 
absorbing gas is not smoothly distributed in space but is collected in discrete clouds with 
small hlling factor, the spectra show occasional narrow features of signihcant optical depth 
even in regions where the mean optical depth is low. The ratio p is useful for predicting 
optical depths averaged over velocity width AV broader than the linewidth of a single cloud, 
< 5 km s“^. 

Using the cool cloud temperature distribution, derived in section 4, together with < k >, 
we can derive an approximate value for the peak brightness temperature of the Galactic 21- 
cm emission. It was noticed as early as the 1950’s in surveys of the Galactic plane that 
the 21-cm emission typically peaks at Tg ~100 to 125 K, regardless of the direction or 
the velocity gradient. This is a curious fact, clearly seen in the spectra of the Weaver and 
Williams (1972) survey, for example. Since we know that the neutral hydrogen is spread 
fairly smoothly from the inner Galaxy to well outside the solar circle, we might expect 
the peak brightness temperature to correspond to low values of the velocity gradient, as 
it naturally would if most lines of sight were optically thin. On the other hand, we know 
that most of the Hi is in the warm medium (roughly 75 %, Kulkarni 1983), so the peak 
brightness temperature does not simply correspond to the cloud temperature, Tcooi, since 
the warm phase gas contributes most of the column density at any velocity. It is clear from 
emission-absorption spectrum pairs at low latitudes, like those in hgures 2-8, that the peak 
brightness in emission does not always correspond to the velocities of high optical depths; 
on the contrary, deep absorption lines often correspond to dips in the emission brightness 
(i.e. HI self-absorption). It is a combination of the opacity, < k, >, the cloud temperature, 
Tcooi, and the warm phase density, Uw, that determines the peak brightness temperature. 

On a line of sight for which p > 1, the 21-cm line is fairly optically thick, and the typical 
line of sight distance. As, needed to achieve average optical depth r = 1 over some velocity 




36 


width AV is 


As 


Ar 

< K> 


( 17 ) 


since 


< K> = 


EW 

As 


Iav ^ dv 
As 


( 18 ) 


If we assume that the warm phase gas has mean density n^, then its column density over 
this line of sight interval is 


AT _ A _ ^10 AV^ /in'! 

< K> 

which, if there were no absorption, gives a velocity integral for the 21-cm line emission of 

[ T^,u dv = ^ (20) 

J 

where the subscript u indicates that this is the unabsorbed value, not what we actually see. 
Assuming once again that this warm phase emission is spread smoothly in velocity due to 
its broad natural line width (> 10 km s“^), we can take 

[ dv ~ AV^ (21) 

Jav 


which gives 


Ch AV^ Ch < fi) > 


( 22 ) 


As in equation 3, we can use Ch = 5.9 x 10“^ to get in K for riyj in cm“^ and < k > 
in km s“^ kpc“^. For a typical value of ?7.^=0.25 cm“^ this gives or about 80 

K if < K >~ 5 km s“^ kpc“^. Note that we have not had to assume a value of AV for this 
calculation, it can be anything between the typical absorption linewidth of <5 km s“^ and 
the typical warm phase linewidth >10 km s“^. 


The brightness temperature we expect from regions with p > 1 is different from this, 
however, since is partially absorbed by the cool clouds before it gets to us, and the cool 
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clouds themselves contribute some emission, given by Tcooi (1 — If the warm phase gas 

is well mixed around and among the cool clouds, then its attenuated brightness temperature 
is given by 


Tw = T^,u - (23) 

T 

Adding on the emission from the cool cloud(s), Tcooi ( 1 — gives the total brightness 

temperature seen in emission 


Tb = Tcooi (1-e-") + (24) 

r 

For our assumed r ~ 1 this gives ~0.6 {Tcooi + Ty^,u)-, or about 0.6 (65K + 80K) = 87K. 

In the alternative geometry assumed in section 4, where the warm phase gas is partly 
in front of the absorbing cloud and partly behind, we get 

Tb = Tcooi (l - e“^) + [l - e (l - e"^)] (25) 

or, in the notation of equations 10 - 14, 

Tb = Tcooi X + T^^u (1 - ex). (26) 

For e=0.5 and a;=0.63 (t=1) this gives 102 K. In both cases the contributions from the warm 
and cool phases are about equal at about 50 K each. Going to higher optical depths does 
not change these numbers much. For r —cx) (x — 1) equation 26 predicts Tb — 115 K and 
in equation 24 the peak brightness actually decreases to approach Tcooi — 65K. Of course 
the actual value of the peak brightness temperature will fluctuate around these values due 
to the relatively wide distribution of values of Tcooi- The highest values we see presumably 
correspond to the warmer clouds {Tcooi ^ 100 K) and regions of relatively low < k, >, i.e. 
long path lengths between clouds. 

In face-on disk galaxies the peak brightness temperature is much less dependent on the 
warm phase gas, since our lines of sight have path lengths of only a few hundred parsecs 
through the neutral hydrogen layer before hitting (or missing) a cloud. This explains why the 
brightness temperatures of the Hi emission from nearby galaxies show such strong variation 
from place to place (Braun 1997). This ’’high brightness network” probably traces the 
structure of the cloud phase of the medium, which may include a significant fraction of the 
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warm gas with it. Payne, Salpeter and Terzian (1983), estimate that only abont 30% of the 
total Hi emission is from the widespread, “interclond” medinm, the rest is associated with 
the clonds, of which abont half must be cool and the rest warm. Braun (1997) sees many 
examples of brightness temperatures above 150 K; the SMC shows a similarly high peak 
brightness temperature of 179 K (Stanimirovic et al., 1999). The implication of the high 
ends of the histograms in hgures 14 and 17 is that there must be some clouds in the Milky 
Way with T^ooi as high as 150 K or more, although they are rare compared with the 50 to 100 
K clouds. If we could get outside the Milky Way and look back, we might see those clouds 
unobscured by the cooler, more common clouds, so that seen from outside the Milky Way 
also may show brightness temperatures of 150 K or more. Braun (1997) also sees a radial 
variation in the brightness temperatures of the “high brightness network”, with brightness 
temperatures increasing with radius. This may be due to a radial variation of Tcooh but it 
could also be explained by a radial variation in < k > similar to what we hnd for the Milky 
Way in section 3; if < k > decreases with Rgai then increases with Rgai in equation 22, 
which makes the peak brightness temperatures increase. 

The recognition that the peak brightness temperature of the 21-cm line at low latitudes 
results from a trade-off between absorption in the cool phase and emission from both warm 
and cool was hrst described by Baker and Burton (1975) in the context of radiative transfer 
models to £t the low latitude survey data then available. The role of the velocity gradient 
in the peak brightness temperature problem is discussed further by Burton (1993). The 
only innovation in the discussion above is to use the measured values of < k > together 
with the mean Hi density to get a rough quantitative estimate for the value of the peak 
brightness temperature, rather than making a full cloudy medium model of the radiative 
transfer like that of Liszt (1983). The opacity function, < K{r) >, incorporates the optical 
depth predictions of these models without making specihc assumptions about clouds sizes 
or geometries. 


5.2. HI Self-Absorption 

Low latitude surveys of the Galactic 21-cm emission like the SGPS and the GGPS have 
turned up many examples of HI self-absorption (HISA, see Gibson et al., 2000, Minter et al. 
2001, and Li and Goldsmith, 2002 for recent reviews). Earlier surveys with large single dish 
telescopes had shown that HISA is relatively common at low latitudes (Baker and Burton, 
1979). It is often, but not always, associated with molecular line emission (Burton et al. 
1978, Knapp 1974). 

HISA is not true self-absorption, but absorption by cool foreground Hi clouds of the 
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background 21-cm line emission from warmer, more distant gas. Although temperature 
determinations for these HISA clouds are imprecise, typically they must be cooler than 40K, 
and sometimes as cold as 10 to 15K to fit the observations. Such cold HI may be inside 
molecular clouds, but this is not always the case. Although we do not yet have a complete 
catalog of such clouds, there are some large examples (Knee and Brunt 2001) which suggest 
that they may be common in the ISM. In this analysis there are some four clouds whose 
apparent temperatures are below 20 K, and several more between 20 and 40 K, depending on 
e. The most obvious is in box 4 toward G326.65+0.59 (-24 < v < -17 km s“^ on figure 13). 
The effect of the cloud is to cause a deep dip in the emission spectrum at these velocities. 

In principle, any cloud whose optical depth profile corresponds to a dip in the emission 
spectrum, rather than a peak, qualifies as HISA. In the notation of equation 14, this criterion 
amounts to 


s = C 2 - e {cq + Cl) <0 (27) 

so that the dependence of Tb on x has a negative coefficient. For the e=0.5 case, the 
median value of s is -2.5; fully 67% of the lines have slope s <0. Thus HISA is in fact more 
common than the alternative, i.e. more than half of all HI clouds, defined as distinguishable 
features in the optical depth spectrum, correspond to dips rather than peaks in the emission 
spectrum. The HISA clouds which can be distinguished on the HI emission survey maps 
must be a relatively small subset of this population, those which are either particularly close 
or particularly cold, so that they stand out in their obscuring effect on the background 
emission. 

At low latitudes and velocities corresponding to the inner galaxy, the HISA effect is 
amplified for clouds on the near side of the sub-central point due to unrelated emission from 
gas at the same velocity on the far side. Depending on the latitude, this should typically 
increase the expected value of e from 0.5 to as high as 0.75. Using e=0.75 changes the 
median value of the slope s above only slightly, from -2.5 to -2.3, but still more than half of 
the clouds (55%) show negative values of s. 


6. Conclusions 

Hi spectra from low latitude surveys of the inner galaxy can only be understood in terms 
of radiative transfer in a partially absorbing medium. The gas is far from optically thin, and 
it is far from isothermal, so neither of these simplihcations can be used to understand the 
brightness temperatures we see in emission. In this paper we study emission-absorption 
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spectrum pairs to try to understand the relationship between the optical depth and the 
brightness temperature through various clouds along several lines of sight. 

This paper concentrates on a few spectra in a relatively small area of the Galactic 
plane in the fourth quadrant. We hope that the techniques described here will be useful to 
analyze spectra toward many more continuum sources throughout the SGPS, and perhaps 
in other similar surveys as well. Since the sample is small, the results described here must be 
considered as only representative; more precise, quantitative measurements of the properties 
of the cool neutral medium will come from larger samples of spectra. 

The three quantities that determine the characteristics of the emission spectra that we 
see in low latitude surveys are the warm phase HI density, r^^„, the temperature of the cool 
clouds, Tcooi, and the average density of the cool medium, which sets the value of < k >. 
This paper discusses the last two, the hrst has been determined from earlier surveys with 
lower resolution (see reviews by Bnrton, 1988, Dickey and Lockman, 1990). With these three 
qnantities in hand, we can understand the peak brightness of the HI emission, the fraction 
of atomic gas in the warm and cool phases, and the signihcance and abundance of HISA 
features. 

The major hndings here are: 

• The mean opacity, < k >, increases as we go inward throngh the Galactic disk, from 
abont 5 km s“^ kpc“^ at the solar circle to as much as 25 km s“^ kpc“^ in the molecular 
ring. This resnlt is in agreement with the finding of Kolpak et ah (2002) for the first 
qnadrant, bnt different from the earlier result of Garwood and Dickey (1989). 

• The temperatnre of the HI in the cool clonds, Tcooi, is lower than what we would 
derive from a simple division of the emission by the absorption, Tgp. This agrees 
qnalitatively (bnt not quantitatively) with results of htting high and intermediate 
latitude emission-absorption spectrnm pairs by Heiles and Troland (2002). The cloud 
temperatures derived here are similar but slightly higher than temperatnres obtained 
for HI clouds in the SMG (Dickey et ah 2000) and in the LMG (Mebold et ah 1997) 
using a similar method. Thongh the median clond temperature is around 65K, the 
distribution of Tcooi has a tail reaching to very cold values (well below 40 K). 

• Putting together the values of < k > and Tcooi we can derive a rongh prediction for the 
peak brightness temperatnre of the Galactic HI emission, that agrees with the ~125K 
that is seen, independent of the velocity gradient, in single dish emission surveys. 

• The slopes of the htted curves in the Tem, (l-e“"^) plane suggest that HISA is a very 
widespread phenomenon. The striking examples of HISA clouds that we notice in 
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emission surveys are only the particularly cold or fortuitously placed cases. 

Extending these techniques to the full survey area will indicate whether these conclusions 
apply to the inner galaxy as a whole, and also give a much better determination of the 
distribution of cool cloud temperatures. 
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A. Least Squares Solution 

The best fit values of Cq, Ci, and C 2 in equation (14) are given by linear combinations of 
moments of the form shown in equation (15). The derivation of these formulae is very similar 
to the problem of finding the coefficients of the best £t polynomial. Begin by defining nine 
combinations of the moments : 


Yi = y - e xy 
Y 2 = vy — e vxy 
Y'i = xy 

Ki = 1 — 2ex + x"^ 

K 2 = V — 2e ¥x + e'^ V x"^ 
Ks = X — e x"^ 

A4 = — 2e X + 

ATs = vx — e V x"^ 

Ae = a? 

then combine these to make four mixed quantities : 


Hi = 



H 2 



Yi K2\ 

Hi ) 
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(AlO) 
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IC 


Ki 


Kl 


H, = [K, - ^]{K, - 


Ko K- 


H, = ^^5 

and finally combining these we get simply 


C2 = 


Hi 
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Ho 


Hi -Ha 


Y2-^ - C2 (K, - 
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K2 Ks 
Kl 


Ka - 


Kl 


and 


Co 


Vi - K2C1 - K3C2 
Kl 
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(AM) 

(A15) 
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B. Comparison with the Method of Heiles and Troland 

Using the Arecibo telescope, Heiles and Troland (2002a, b) have recently carried out a monu¬ 
mental survey of emission and absorption in the 21-cm line toward some 75 bright continuum 
sources at high and intermediate Galactic latitudes. This survey is complementary to the 
low latitude surveys because it samples primarily gas within two kpc of the sun, but it sam¬ 
ples this nearby gas very well, since the continuum sources are bright and even with the 
3' (FWHM) beam of the Arecibo dish the spatial resolution at the distance of the gas is 
typically less than one parsec. We have analysed their spectra using the method described 
in section 4.3 above, to hnd the distribution of cool phase temperatures shown on hgure 18. 
Counting each cloud (absorption feature) equally gives median temperature 31 K for this 
distribution. Weighting by the number of velocity channels in each cloud raises this median 
temperature to 35 K. 

This result is somewhat cooler than the median temperature found by Heiles and Troland 
(48K), so it is worth looking in detail at the difference between their analysis and ours in a few 
pairs of spectra. Table 4 from Heiles and Troland (2002a) gives a comparison between their 
method and a simpler version of our technique, which they call the “slope method” (Mebold 
et ah 1997). For two sources they take e ~ 0.5 on this table, 3C333 and PKS0428-I-20 
(which they correctly identify by this name on their table 2). For 3C333 they measure 29 
K for Tcooi from the slope method, and 25 K by htting Gaussians, whereas we measure 26 
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K using method 2 of section 4.3 above. For PKS0428+20 they measure 41 K and 2 K from 
the slope method for the two dominant features centered at +3 and +10 km s“^ respectively 
(hgures 19 and 20), while we hnd 38 K and 27 K. From Gaussian htting, Heiles and Troland 
hnd 29 K for the hrst and two blended components at 27 K and 34 K for the second. These 
differences are minor, and they probably result from the apportionment of the warm phase 
emission between the warm and cool components of Heiles and Troland’s Gaussian analysis. 

The analysis method used by Heiles and Troland is to hrst decompose the absorption spec¬ 
trum into the sum of several Gaussians (three in the case of PKS0428+20), and then to 
decompose the emission spectrum into the sum of several more Gaussians (hve in this case), 
of which the centers and widths of the hrst three are set by the values obtained from htting 
the absorption. This technique accomplishes the all-important separation of the emission 
into warm and cool components, where the sum of the emission Gaussians that are not seen 
in the absorption correspond to in equation 13. 

There is no ’’right answer” to the problem of htting multiple Gaussians, the result depends 
on the number of components and the hrst guesses for their parameter values. The linear 
least squares method has this advantage, that the solution for Tcooi and is unique given 
only an assumed value of e and the velocity ranges to use for the htting. At low latitudes, 
the approach of Heiles and Troland is not an option, since the emission spectra are so 
complicated that Gaussian htting is hopeless. In the higher latitude directions studied in 
the Arecibo survey it is possible to decompose the emission as well as the absorption into 
a few Gaussians, and this gives more physical information about the gas than the simple 
two-phase approximation used in this paper. In general the contribution to the emission by 
the components that are seen in absorption is less than the contribution from those that are 
not seen in absorption, i.e. the WNM components; typically more than half of the emission 
integral (column density) comes from Gaussians that are not visible in absorption, and so 
must be warmer than 150 to 200 K. There is a good qualitative agreement on this result 
between this paper and those of Heiles and Troland, in spite of some minor disagreements 
on specific values of Tcooi given by the two methods. 

All methods of analysis seem to agree on the result that the typical temperatures of absorbing 
clouds in the high and intermediate latitude sample observed from Arecibo are cooler than the 
inner galaxy clouds from the SGPS studied in this paper. This may reflect a real gradient in 
the cool phase temperature, either with Galactic radius, R, or with height above the plane, 
z. This may result from the different conditions inside cloud complexes, i.e. the larger 
environments containing cool atomic and molecular clouds and HI I regions. The Arecibo 
sample contains HI clouds that are mostly not associated with H H region complexes, they 
are in a typical interstellar radiation field. If more of the Hi in the inner galaxy is recent 
photodissociation product due to the greater abundance of hot stars as suggested by Allen 




Fig. 18.— Distribution of cool cloud temperatures, Tcooi, for the absorption lines in the 
survey of Heiles and Troland (2002a). The median cool phase temperatnre is 31 K. 
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Fig. 19.— Details of the htting of the emission toward the lower velocity absorption compo¬ 
nent of PKS0428-I-20. Heiles and Troland find 29 K for the cool phase temperatnre, while 
we hnd 38 K. 
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(2001), then the typical cool phase temperature might be elevated, or there might be still 
more blending of warm, cool, and intermediate temperatures that could bias our results 
to higher values of Tcooi- If further study of the SGPS data conhrms that there is a trend 
toward lower cool phase temperatures in HI clouds with increasing R, then the corresponding 
decrease in < k > with R suggests that the abundance of cool clouds is decreasing rapidly 
with increasing R. 
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Fig. 20.— Details of the fitting of the emission toward the higher velocity absorption com¬ 
ponent of PKS0428-I-20. Heiles and Troland decompose this into two Gaussians with tem¬ 
peratures 27 K and 34 K, while we hnd 27 K for the combination. 
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Table 2. Fit Parameters, e = 0.5 


source 

component 


Vo,i 

km s“^ 

cr* 

km s“^ 

Vi 

km s“^ 

Vi+l 
km s“^ 

Co 

K 
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Cl 

K 

C2 

K 

14 

C2 

K 

T^o{vi) T^{vi+i) 

K K 

method 3 

km s~^ 

Lod 2, eq. 

G331.35+1.07 

1 

1.65 

-82.1 

3.6 









G331.35+1.07 

2 

0.33 

-71.4 

4.2 









G331.35+1.07 

3 

0.55 

-54.9 

6.5 









G331.35+1.07 

4 

0.81 

-30.4 

2.0 









G331.35+1.07 

5 

0.52 

-24.5 

21.3 









G331.35+1.07 

6 

0.54 

-5.8 

3.8 








1 

G331.35+1.07 

7 

1.00 

4.3 

3.1 









G331.35+1.07 

8 

0.19 

18.8 

2.9 








1 

G331.35+1.07 

9 

0.32 

37.5 

3.5 









G326.45+0.90 

1 

1.49 

-43.3 

3.8 

-45.8 

-36.7 

25.33 

65.90 

63.24 

80.83 

-1.43 

18.34 

G326.45+0.90 

2 

0.87 

-40.3 

1.5 

-45.8 

-36.7 

25.33 

65.90 

63.24 

80.83 

-1.43 

18.34 
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3 

1.26 

-31.3 

2.0 

-33.4 

-28.5 

33.37 

61.77 

62.20 

60.45 

0.10 

34.58 
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4 

1.12 

-17.2 

3.0 

-20.3 

-13.7 
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37.51 

1.23 
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5 

0.32 

-4.5 

3.8 

-7.9 

-2.1 
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53.51 
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1 
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-3.65 
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1.32 
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3.7 
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4.99 

G326.65+0.59 
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1.46 
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-5.4 
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13.18 
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0.26 
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3 
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-65.6 
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107.65 
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107.4 
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4 

1.69 

-46.7 

12.1 

-54.1 

-37.6 

36.99 
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5 

1.33 

-37.8 

2.7 

-37.6 

-34.3 

40.17 
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65.23 

79.1( 
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6 

0.83 

-16.0 

4.6 

-22.7 

-7.9 

46.11 

46.66 

42.22 

43.6! 

G328.31+0.45 

7 

0.40 

-1.5 

6.8 

-9.5 

0.4 

37.67 

42.22 

70.10 

39.0i 

G328.42+0.22 

1 

0.08 

-118.2 

1.6 

-119.2 

-115.1 

2.44 

4.21 

9.57 


G328.42+0.22 

2 

0.52 

-107.5 

2.4 

-109.3 

-104.3 

26.65 

34.94 

57.20 


G328.42+0.22 

3 

0.99 

-101.6 

2.7 

-102.7 

-99.4 

28.87 
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4 

2.17 

-87.8 

13.3 

-96.9 

-71.4 
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5 

1.38 
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54.2! 
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6 

1.60 

-51.0 

9.2 
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G328.42+0.22 

7 

1.60 
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114.9 
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-31.8 
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3 
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4 
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5 
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1 
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2 
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3 
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4 
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5 
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6 
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7 
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8 
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1 
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T),i ^0,i 

km s“^ km s“^ 


0.82 

-17.7 

5.3 

1.60 

-5.2 

8.6 

1.60 

-5.2 

8.6 

0.51 

27.7 

3.8 

1.75 

-89.3 

4.7 

1.80 

-63.8 

7.1 

2.23 

-47.1 

6.8 

0.33 

-22.8 

7.9 

0.52 

-0.6 

4.3 

1.59 

-87.1 

3.2 

0.32 

-76.4 

2.8 

1.36 

-65.6 

2.5 

1.47 

-60.3 

2.4 

11.95 

-50.4 

2.3 

1.85 

-45.4 
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0.56 

-31.8 

1.8 

0.37 

-17.0 

9.8 

0.98 

-1.6 
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3.96 

-91.7 
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1.68 
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4.0 
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112.66 

17.03 


23.4 

31.7 

24.08 

19.41 

35.30 

15.81 



-99.4 

-70.5 

-54.1 


-83.7 

-64.8 

-45.0 


15.30 72.91 43.43 

22.64 69.36 42.01 

3.33 38.17 70.57 


103.9 
-8.4( 

107.9 














Table 2—Continued 


source 

component 


V0,i 

km 

km 

G331.52-0.07 

4 

2.90 

-53.4 

12.0 

G331.52-0.07 

5 

0.75 

-18.3 

2.3 

G331.52-0.07 

6 

0.27 

-1.6 

7.6 

G331.52-0.07 

7 

0.27 

-1.6 

7.6 


Vi 

Vi+l 

Co 

Cl 

C2 

C2 

km 

km 

K 

K 

km s~^ 

K 

K 



method 2, eq. 

14 


-44.2 

-32.6 

6.05 

70.57 

53.99 

114.01 

-21.9 

-12.8 

11.02 

44.27 

36.73 

48.54 

-11.2 

-2.9 

-17.04 

36.73 

75.41 

33.89 

-2.1 

8.6 

-3.15 

75.41 

23.57 

75.65 












